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ABSTRACT 



Context. Combined spectroscopic abundance analyses of stable and radioactive elements can be applied for deriving stellar ages. The 
achievable precision depends on factors related to spectroscopy, nucleosynthesis, and chemical evolution. 
Aims. We quantify the uncertainties arising from the spectroscopic analysis, and compare these to the other error sources. 
Methods. We derive formulae for the age uncertainties arising from the spectroscopic abundance analysis, and apply them to spectro- 
scopic and nucleosynthetic data compiled from the literature for the Sun and metal-poor stars. 

Results. We obtained ready-to-use analytic formulae of the age uncertainty for the cases of stable+unstable and unstable+unstable 
chronometer pairs, and discuss the optimal relation between to-be-measured age and mean lifetime of a radioactive species. 
Application to the literature data indicates that, for a single star, the achievable spectroscopic accuracy is limited to about ±20 % 
for the foreseeable future. At present, theoretical uncertainties in nucleosynthesis and chemical evolution models form the precision 
bottleneck. For stellar clusters, isochrone fitting provides a higher accuracy than radioactive dating, but radioactive dating becomes 
competitive when applied to many cluster members simultaneously, reducing the statistical errors by a factor V5v. 
Conclusions. Spectroscopy-based radioactive stellar dating would benefit from improvements in the theoretical understanding of 
nucleosynthesis and chemical evolution. Its application to clusters can provide strong constraints for nucleosynthetic models. 

Key words. Sun: abundances - Stars: abundances - Stars: evolution - Stars: fundamental parameters 
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The determination of absolute ages for stars is of paramount 
importance in astrophysics since ages allow the timescales on 
which stellar populations and their chemical composition evolve 
to be established. The presence of long-lived unstable elements 
in stellar atmospheres offers a way to date stars by investigating 
their radioactive decay spectroscopically. In its simplest form, 
one determines the photospheric abundance of a pair of ele- 
ments, of which one is unstable. By means of theoretical mod- 
eling of the nucleosynthetic production channels for the two el- 
ements, the original abundance of the unstable nucleus can be 
estimated from the abundance of the stable one. This allows the 
time-span over which the unstable element decayed to be deter- 
mined. 

The idea emerged early after the discovery of the rapid neu- 
tron capture process (short r-process) as the source of produc- 
tion of radionuclides such as thorium and uranium (see Cowan 
et all 1991, and references therein for an account of its early ap- 
plications). Since then it has been applied i n a number of cases 
(among the most recent Fre bel et al.ll2007t) but, typically, as a 
"side product" of analyses of the more general n-capture enrich- 
ment patterns. To our knowledge, the only effort in recent years 
to exploit radioactive dating to systematically d erive the age of a 
signific ant sample of stars has been the one bv ldel Peloso et alTI 
(l2005allbllcr) . 

One interest in such dating techniques lies in its applicabil- 
ity to field stars for which no distance information, hence no ac- 
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curate information on surface gravity or absolute luminosity, is 
available, so that they cannot be placed with confidence on the- 
oretical isochrones in the Hertzsprung-Russell diagram (HRD). 
For cluster stellar populations, on the other hand, HRD fitting 
usually provides rather precise absolute ages (e.g. Vandenberg 
et al. ll996l) . However, for stellar clusters (especially globular 
clusters), radioactive dating could produce an age estimate of 
comparable or even higher precision if applied to a large enough 
sample of stars, and would as such be very valuable as a test of 
the theoretical stellar isochrones. 

Although potentially effective, spectroscopy-based radioac- 
tive dating is difficult to apply. The most easily measured long- 
lived radioactive element, thorium, is only accessible (with the 
exception of strongly r-process enriched stars) through a single, 
weak line in the optical spectrum of stars, in which high accuracy 
is needed in the measurement to provide a reliable age. Besides 
lying on the red wing of a Fe-Ni blend, this unique Th n line is 
itself also blended by weaker lines of vanadium and cobalt. The 
relative st rength of the variou s components for the Sun are dis- 
cussed in Caffau et al. (2008) where the line-list of del Peloso 
et al. (l2005al) was applied. Although rather easily measured, the 
most frequently used reference element, europium, is believed to 
constitute a less-than-optimal reference. In fact, its nucleosyn- 
thetic ties with Th are rather weak, and the production ratio is 
sensitive to the originating supernova characteristics (Rratz et 
al. l2007tlFarouqietal.ll2008h . 

Besides the astrophysical problem that the r-process site has 
not been unambiguously identified, it is unclear how universal 
the resulting abundance pattern of the r-process nucleosynthesis 
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actually is. In view of the complex nuclear reactions shaping the 
r-process it is not obvious that the outcome is the same for a po- 
tentially wide range of the physical conditions governing the for- 
mation site(s") dGorielv & Clerbaux Hl999l : lArnould et al. II2007L 
ISneden et al.l u008: Roe derer et al.ll2009h . To minimize the influ- 
ence of poorly constrained physical conditions at the formation 
site it is desirable to use a reference element with an atomic mass 
as similar as possible to thorium. This should reduce the uncer- 
tainties in the formation properties to the uncertainties related to 
the properties of the involved nuclei as such, and environmental 
factors sho uld cancel out. As pointed o ut before by various au- 
thors (e.g.. iGorielv & Clerbaux I U999) the Th/U pair is in this 
sense optimal, however, again difficult to app ly since uranium is 
difficult to detect spectroscopically. Recently. iRratz et ail d2007l ) 
suggested the element hafnium as suitable reference element for 
thorium-based radioactive dating. Hafnium is closer in mass to 
thorium than europium, but more importantly, Kratz and collab- 
orators showed by detailed modeling that the production ratio 
Th/Hf is not strongly varying over a wide range of neutron fluxes 
during the r-process. Since hafnium is a spectroscopically acces- 
sible element this makes it an interesting candidate as reference 
element. 

Observationally, the forthcoming introduction of Extremely 
Large Telescopes (ELTs) is likely to make radioactive dating of 
stars a much more exploitable technique, since the high signal- 
to-noise (S/N) ratio required is one of the main difficulties asso- 
ciated with such studies. In this perspective, a precise assessment 
of the uncertainties associated with such technique is currently 
lacking, and attempted here. We derive the formulae for the pre- 
cision to which spectroscopy can provide age estimates. They 
are subsequently applied to abundance analyses of stars com- 
piled from the liter ature. This paper is a "spin-off" of our work 
(ICaffauetal.12 008) on the thorium and hafnium abundance in the 
Sun, and the Th/Hf chronometer pair often serves as an example 
in the following sections. 

2. Formulae related to radioactive dating 

2.1. Age uncertainties 

Any sample of radioactive nuclei decays according to the expo- 
nential law: 



n x (t) = nxoexp(-f/r) 



(1) 



where n X o is the number of nuclei at ?=0, and nx the number of 
remaining nuclei at time t; r is the mean lifetime. Explicitly, t is 
given by 



t = rln 



nxo 
n x 



(2) 



Knowing the initial and present number of nuclei, one can obtain 
the time interval that has elapsed since t =0. In practice, to use the 
decay processes to determine an age, one measures the number 
of nuclei of the unstable species (nx) relative to the number of 
nuclei of a stable element which is related to nxo- This is, e.g., 
the case for the Th/Hf pair. Let us call n« the number of nuclei of 
the stable reference species related to nxo. Let us further assume 
that the two nuclei have a common nucleosynthetic origin and 
that the relation between the initial number of nuclei nxo and n R 
is given by 



nxo = fin R 



(3) 



where /3 is the production ratio between the respective numbers 
of nuclei. The production ratio may be predicted by theoretical 



considerations (see Kra tz et ai1l2007L and references therein). 
From Eq. (0 and Eq. (0 one can write 



rln 



nx 



(4) 



All astrophysically relevant unstable elements are trace 
species. Let us therefore suppose that the lines corresponding 
to the transitions we are interested in, are weak. In this case, the 
relation between equivalent width (W hereafter) and abundance 
(n) is linear: 



W = C n or log W = log C + log «, 



(5) 



C being a constant defined by the line transition under consider- 
ation (dependent on oscillator strength and excitation potential) 
and the thermal structure of the stellar atmosphere under inves- 
tigation. 

Many lines - including the Th n resonance line - are blended 
with other species, and in practice the abundance is obtained by 
fitting a synthetic profile. In this case the equivalent width is ob- 
tained theoretically from the best-fitting profile and the underly- 
ing abundances. The error on the equivalent width may then be 
obtained by propagating the error in Eq. 0. One may also try 
to derive an equivalent width for the blended line by subtract- 
ing the theoretical equivalent width of the blending components 
from the total W of the blend. However, in general, this ignores 
saturation effects, and thus the "corrected" W differs from the 
true one. As an example we may consider the Th n resonance line 
in the Sun. Let us consider the theoretical computation of Caffau 
et al. (120081) based on a 3D CQ 5 BOLD model: the W of the 
blend Th n+Co i+V n, for an assumed abundance A(Th)=0.09 
is 0.74pm, the W of just the Coi+Vii blend is 0.31pm, while 
W(Th) = 0.46 pm. In spite of the fact that all the lines involved 
are very weak, it is clear that the "corrected" W (0.43 pm) differs 
from W(Th), due to saturation effects. Thus one should always 
use the equivalent width of the line, rather than the "corrected" 
W, although in many cases the difference is so small that it may 
be safely ignored. Alternatively, one may work with Eq. (0 and 
use the obtained abundances directly. In this case one has to es- 
tablish uncertainties by Monte-Carlo simulations. 

Substituting Eqs. (0 and (0 into Eq. (0 we obtain 



rlnL0 



Wr Cx 
W x C R 



= t \ ln/3 + In W R - In W x + In | ^ 



(6) 



The achievable precision in t is hence limited by the uncertainty 
in the knowledge of /?, errors in the measurement of Wx and 
Wr, as well as by systematic errors in modeling the stellar at- 
mosphere and line transitions, i.e. errors of the constants Cx and 
C R . 

From the standard formula for the propagation of uncorre- 
lated errors we get for the variance of t 



O-(t) = T < 



ajW R ) 

Wr 



ajWx) 
W x 



cr(C x /C R ) ' 
Cx/C R / 



1/2 



(7) 



This equation is valid under the presumption that the errors of /?, 
W x , W R , and (Cx/C R ) are uncorrected. The errors of Cx and C R 
themselves may often be more or less correlated. In fact, Th and 
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Hf abundances are derived from singly ionized species which 
constitute the dominant ionization stage in the late-type stars of 
interest here. Moreover, the diagnostic lines of Th and Hf emerge 
from low lying energy levels which makes the Boltzmann factor 
and consequently the temperature sensitivity of the level popula- 
tions modest. The similar atomic properties of Th and Hf lead to 
similar formation regions in the photosphere. Hence, all system- 
atic errors due to uncertainties in gravity, effective temperature 
and even the details of the thermal structure of the employed 
model atmosphere should largely cancel out. Hf is the heaviest 
element prior to Th whose observable spectroscopic signatures 
are this favorable. What is left is the systematic error related to 
the ^/-values. 

Taking the Sun and the Th/Hf pair as an example, we have 
o-(W R )/W R ~ a 0.03/ 0.3 = 0.1 (see Hf n A 409.3 nm in Table 3 in 



ICaffau et alJ (l2008h\ and gW x)/W x ~ 0.03/0.4 = 0.075 (see 



Table 4 in lCaffau et alj |2008). Since 232 Th has a mean lifetime 
of t=20.3 Gyr, we obtain from Eq. (0 under the most optimistic 
conditions, assuming cr(J3) = and ct(Cx/Cr) = 0, an error in 
the age determination of the order of cr(f)~2.5 Gyr. 

Even this rather limited precision is probably beyond the 
present state-of-the-art in stellar spectroscopy. However it is not 
implausible to reach a precision of the order of 2 Gyr in the fu- 
ture. The uncertainties due to the S/N in the spectrum may be 
considerably reduced by the next generation of ELTs, and uncer- 
tainties in the blending lines can be reduced through laboratory 
measurements and/or theoretical calculations. The precision cur- 
rently attainable on the dating of Globular Clusters, using Main 
Sequence fitting, is a factor of two better (about 1 Gyr), how- 
ever it may not be applied to field stars. For field stars which are 
evolved off the Main Sequence, if the distance is known, an age 
can be derived through comparison with theoretical isochrones. 
However, distances are difficult to measure and a small error in 
the distance propagates to a great error in age. In principle, ra- 
dioactive dating has the potential to provide ages of individual 
field stars, even if their distances are not precisely known. 

Unfortunately, in the literature equivalent widths and their 
uncertainties are usually not given so that a direct application of 
the previous formulae is not possible. In Sect. l2.2l we will see that 
they are nevertheless useful to derive a statement of the optimal 
decay time of a species for measuring a particular time interval. 
For now, we proceed by listing the previous relations in terms 
of the logarithmic abundances A(X) = log ex = log«x/«H + 12. 
Equation (0 for the age becomes 



t = rlnlO {log/? + A(R) - A(X)} , 
with the corresponding uncertainty 

a(t) = rlnlO {<x 2 [log/3] + cr 2 [A(X)] + cr 2 [A(R)] 

- 2a[A(X)] a[A(R)] Cor [A(X),A(R)] }' /2 



(8) 



(9) 



Cor [., .] denotes the linear correlation coefficient. In deriving 
Eq. (0, it was assumed that /3 is uncorrected with A(X) and 
A(R). Relation (0 is consistent with Eq. as long as the statis- 
tical bias due to the change from linear to logarithmic quantities 
is small, i.e. that for any involved quantity q 



o-(q) 
q 



cr(ln q) = In 10 cr(log q). 



(10) 



Equation (0 now explicitly contains the correlation among the 
measurements of the abundances of the two involved species. If 
abundance errors are positively correlated - as discussed previ- 
ously and the Th/Hf pair serving as example - this can lead to a 



substantial reduction of the overall error with respect to the un- 
correlated case. However, care has to be exercised when apply- 
ing the formulae to literature data (as we intend to do). Quoted 
abundance errors are often purely statistical, uncorrected un- 
certainties while correlated uncertainties related to atmospheric 
models are specified separately by providing sensitivities against 
changes of the stellar parameters. 

2.2. Optimal decay times 

We now discuss how the achievable precision depends on the 
time interval t we want to measure and the decay time of the ra- 
dioactive species, t. From Eq. or (0 it might appear that a 
radioactive species with short mean lifetime r provides the high- 
est precision on t. However, this reasoning ignores the fact that a 
rapid decay (r <K t) makes it more difficult to measure Wx pre- 
cisely. To make this connection manifest we eliminate W x and 
rewrite Eq. as 



t 



T ^cr(U A .) 



w R 

2 



1 + 



Cr 



g(C x /C R ) 
CxlC R 



PCx 

2^/2 



exp(2f/r) 



where we have used the relations 

Wr Cr 

W x PC X 
and 

a(W H )=a(W x ). 



(11) 



(12) 



(13) 



Eq. (112b is obtained by combining Eqs. 0, (0, and (0, while 
Eq. ( 113b is a consequence of Cayrel's formula dCavrell [l988) 
which, for a spectrum that is photon noise dominated, states 
that the uncertainty of the W measurement cr(W) depends on the 
FWHM of the line, the pixel size, and the S/N of the spectrum. 
For spectra of moderate resolution (R< 80000), the FWHM is 
essentially the same for all lines since they are not fully resolved 
and the FWHM is simply set by the spectrograph's resolution. 
This is also true for higher resolution spectra if the atoms under 
consideration do not have a very large difference in mass, and 
therefore thermal broadening, or if the dominant broadening is 
non-thermal. In the case of the Th/Hf pair, the mass difference 
is not large, and in late-type stars the macro-turbulent plus ro- 
tational broadening is often of the same order of magnitude or 
greater than the thermal broadening. We may therefore assume 
that <t(W) is the same for any line in a given spectrum. Assuming 
again cr(fi)=Q, and ct(Cx/Cr)=0, and noting that ct(Wr)/Wr is 
time independent, the relative error of t becomes 



£t0 
t 



<T(W R ) t 

t 



W, 



R 



0iW R ) It \ Cr 

F - ; p\ with p = . 

Wr 1 PC X 



l + exp(2,/ T ) 

Cr 



1/2 



(14) 



For brevity, we introduced an auxiliary function F of argument 
and parameter p in Eq. (TBI) . The smaller F at given p, the 
higher the precision of the age determination. The best (relative) 
accuracy is obtained for the time where the function F shows 
a minimum. As shown in Fig.0 the function has a minimum 
near t » r. Left of the minimum (f < t), F decreases slowly 
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Fig. 1. The function F^;pj defined in Eq. ( 1141 ) plotted as a 
function of - for different parameters p (solid curves). The 
dashed line shows the case of the Th/Hf chronometer pair un- 
der solar conditions (Thn A 401.8 nm and Hf n A 409.3 nm line 
pair, p=0.26). The the dotted and dash-dotted vertical lines mark 
two characteristic ages, 5 and 12 Gyr, taking the Th (20.3 Gyr) 
and U (6.5 Gyr) mean lifetimes as reference. 

with time, while it increases more steeply beyond the minimum 
(f > r). We can conclude that an optimal relative precision for 
radioactive dating is obtained when one measures a time com- 
parable to the mean lifetime of the unstable element - a perhaps 
intuitively expected result, put on a quantitative footing here. Of 
course, our analysis is idealized by assuming a single, isolated 
line for both the stable and decaying species. In practice, often 
several lines are considered and line strengths are measured via 
spectral synthesis in order to account for line blending. However, 
we think that the overall picture in this more complex situation 
does not change qualitatively. 

We finally consider radioactive dating when measuring two 
lines from two transitions of two unstable elements which we 
call X and Y. There are two decay laws to be considered for the 
two elements: 

«x = n X oexp(-t/T X ) and n Y = n Y0 exp (-t/ry) . (15) 

We suppose that tx < Ty, and that the two elements are, at least 
partially, formed in the same nucleosynthetic process, so that the 
abundances of nuclei of the two elements at the initial instant of 
time are related by nxo = y6«ro- In this way, we obtain 

n x =y8nyexp(-f/-r), (16) 

once t is defined as 

I -1-1. (17, 

T T X Ty 

We then find the same result for cr(t)/t as obtained above in 
Eq. (TBI) , except that r now depends on the mean lifetimes of 
the two elements considered according to Eq. (117b . and the ref- 
erence element R is identified with the longer lived species Y at 
its present abundance. 

In practice, the only two radioactive nuclei which have been 
used to date stars are 232 Th and 238 U with mean lifetimes of 
20.3 Gyr and 6.5 Gyr, respectively. From Eq. ( TTTb the effective r 
for the Th/U pair amounts to 9.6 Gyr. From Fig. [JJit is clear that 
the Th/U pair as well as U/R pairs (R for any stable r-process el- 
ement) are well suited for measuring great ages around 12 Gyr, 



with Th/R being competitive depending on the actual value of 
the p parameter. Towards ages of a few Gyr the sensitivity of 
U/R pairs is about 2 to 3 times higher than the one of Th/R pairs, 
again depending on p. We repeat that only the decay times enter 
the considerations here; it is usually more difficult to measure 
U than Th abundances so that the advantage of the shorter life 
time of U is offset. Finally, we ask whether there are other suit- 
able naturally occurring radioactive elements with lifetimes even 
more favorable for measuring ages in the range of a few Gyr. The 
answer is 'no' : the lifetimes are either much too short or much 
too long - besides further complications like exceedingly low 
abundances of the respective isotopes or the availability of suit- 
able spectral lines. 

3. Ages of metal-poor stars 

As an application of the formulae we derived for the errors on 
age, we consider metal-poor stars for which the thorium abun- 
dance has been measured. For all these stars the abundance of 
at least one other r-process element is available (see Table lA.fi ). 
Using various production ratios /3 available in the literature (see 
Table IA.2I ). we computed the ages of these stars using Eq. ([8} 
and their uncertainties with Eq. ((9). Correlations were neglected; 
correlations due to uncertainties in stellar parameters are usu- 
ally not included in the error budget, and correlations emerging 
from the data reduction (e.g. background reduction) not speci- 
fied by the authors. Our results are summarized in Table [JJ For 
star HE 1523-0901 t he abundance s of t he single elements were 
not available, but in iFrebel et al.l d2007h abundance ratios are 
reported, and measurement errors on the individual abundances 
which is sufficient for our purposes. By using different produc- 
tion ratios we intended to separate the systematic uncertainties 
related to nucleosynthesis and chemical evolution from the spec- 
troscopic uncertainties as such. Since the production ratios are 
treated separately here their related uncertainty in the evaluation 
of <x(f) with Eq. (|9|) is set to zero, cr 2 [log/?] = 0. 
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Table 1. Ages and age uncertainties cr A6C of metal-poor stars and the Sun obtained for different chronometer pairs X/Y. 



X/Y Object 


Age a 


0~ Age 




[Gyr] 


[Gyr] 



Th/Eu 


HD 115444 


11.7 


10.4 


12.3 


14.6 


6.9 


Th/Os 


HD 115444 


30.1 


24.3 






7.7 


Th/Eu 


CS 22892-052 b 


14.8 


13.4 


15.4 


17.6 


6.4 


Th/Hf 


CS 22892-052 


27.7 








7.4 


Th/Os 


CS 22892-052 


19.6 


14.8 






18.5 


Th/Eu 


CS 22892-052 


15.4 


13.1 


15.0 


17.3 


5.3 


Th/Os 


CS 22892-052 


24.4 


18.6 






5.1 


Th/Ir 


CS 22892-052 


25.8 


19.1 


25.2 




12.9 


Th/Eu 


CS 22892-052 


12.1 


11.7 


13.7 


15.9 


4.8 


Th/Hf 


CS 22892-052 


24.6 








6.1 


Th/Os 


CS 22892-052 


26.0 


20.3 






7.7 


Th/Ir 


CS 22892-052 


19.0 


13.3 


19.4 




6.1 


Th/Eu 


BD +17 3248 


7.6 


6.3 


8.3 


10.5 


5.3 


Th/Os 


BD +17 3248 


27.7 


22.0 






6.1 


Th/Ir 


BD +17 3248 


20.7 


14.0 


19.1 




7.7 


Th/Eu 


HD 221170 


11.7 


10.4 


12.3 


14.6 


4.2 


Th/Hf 


HD 221170 


21.9 








2.9 


Th/Os 


HD 221170 


27.1 


22.3 






5.3 


Th/Ir 


HD 221170 


20.7 


14.0 


19.1 




6.1 


Th/Eu 


CS 31082-001 


0.9 


-1.5 


0.5 


2.7 


2.9 


Th/Hf 


CS 31082-001 


15.2 








8.0 


Th/Os 


CS 31082-001 


17.9 


12.1 






5.9 


Th/Ir 


CS 31082-001 


6.5 


-0.2 


5.9 




4.1 


U/Th 


CS 31082-001 


16.6 


13.3 


14.6 




5.9 


Th/Eu 


HE 1523-0901 


11.4 


9.5 


11.0 


13.2 


3.1 


Th/Os 


HE 1523-0901 


16.9 


10.1 






3.1 


Th/Ir 


HE 1523-0901 


18.0 


12.3 


17.4 




3.1 


U/Th 


HE 1523-0901 


15.0 








2.6 


Th/Eu 


M51 K341 


11.0 


9.7 


12.6 


13.9 


6.9 


Th/Eu 


M51 K462 


14.1 


12.7 


14.7 


16.9 


6.9 


Th/Eu 


M51 K583 


6.6 


4.3 


7.2 


8.5 


6.9 


Th/Hf 


M4L1411 


32.5 








3.4 


Th/Hf 


M4L1411 


32.5 








3.4 


Th/Hf 


M4L1501 


25.5 








3.4 


Th/Hf 


M4L1514 


37.2 








3.4 


Th/Hf 


M4 L2406 


32.5 








3.4 


Th/Hf 


M4 L2617 


23.2 








3.4 


Th/Hf 


M4 L3209 


30.2 








3.4 


Th/Hf 


M4L3413 


23.2 








3.4 


Th/Hf 


M4L4511 


27.8 








3.4 


Th/Hf 


M4L4611 


37.2 








3.4 


Th/Hf 


M5 IV-81 


18.5 








3.4 


Th/Hf 


M5 IV-82 


18.5 








3.4 


Th/Eu 


HD 108577 


7.8 


6.2 


8.4 


10.1 


3.4 


Th/Eu 


HD 115444 


9.6 


8.1 


10.3 


11.9 


3.6 


Th/Eu 


HD 186478 


16.6 


15.1 


17.3 


18.9 


4.3 


Th/Eu 


BD +82856 


7.3 


5.8 


8.0 


9.6 


5.5 


Th/Eu 


M92 VII- 18 


6.8 


5.3 


7.5 


9.1 


5.3 


Th/Eu 


Sun 


3.3 


1.7 


3.9 


5.6 


2.0 


Th/Hf 


Sun 


22.3 








2.3 


Th/Os 


Sun 


14.1 


8.6 






4.9 


Th/Ir 


Sun 


10.1 


4.0 


9.6 




4.9 



a Ages and age uncertainties are obtained from formulae f5] and 191 using the abundances listed in TablelA.ll The four "Age" columns give ages based on production rations Q of Kratz et 
al. |2007|) . |Sneden et a l. \ 2003). Schatz et al. 1 20021), and lCowan et alj (20021) (in that order) according to Table |A.2| The Solar a ges are calculated applying r-process fractions as given in 
TableHPl 

b Three independent measurements were obtained for this star which is the reason why certain chronometer pair/object combinations are listed multiple times. 
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We like to point out that we applied the spectroscopy-related 
uncertainties of the abundances as given by the authors. There 
is no general consensus about the appropriate way to estimate 
the error on abundances, and the approach varies from author to 
author. For those elements which have several measurable lines 
the error can be estimated by looking at the RMS line-to-line 
abundance scatter. One further might divide this error by V^V 
where N is the number of lines, provided the errors on each line 
are independent. However, this assumption is often questionable, 
since in most cases the lines are measured from the same spec- 
trum. Any error which affects the whole spectrum (e.g. an error 
in the background subtraction) will affect all the lines in the same 
(or at least in a similar) way. We thus argue that it is more real- 
istic and conservative not to invoke this factor 1 / y/N. 

The ages of the individual metal-poor stars in Table [T] show 
a great spread when considering different r-production ratios, 
and chronometer pairs. Also the age uncertainties stemming 
from the spectroscopically determined abundances, <x Age , are 
substantial. We note that our spectroscopy-related uncertainties 
are commonly greater than the estimates given by the authors 
themselves. Figures [2] and [3] provide a graphical representa- 
tion of the situation for the stars CS 22892-052 and CS 31082- 
001, respectively. We picked these objects since spectroscopic 
measurements for a greater number of chronometer pairs ex- 
ist. CS 31082-001 is of particular interest since it is strongly 
r-process element enhanced allowing the measurement of the 
Th abundance exceptionally from several lines making it more 
reliable. For CS 22892-052 we left out the age estimates with 
the highest uncertainties, stemming from Th/Os and Th/Ir, since 
they were later improved substantially. Both figures illustrate 
that the dispersion due to different production ratios is often 
greater than the dispersion due to abundance uncertainties. This 
is even the case if one restricts the comparison to production ra- 
tios originating from the same author. As already pointed out in 
th e original paper o n the abundance analysis of CS 31082-001 
bv lHill et al.1 d2002h the age derived from the Th/Eu pair is ob- 
viously unrealistic despite rather accurate abundance measure- 
ments. The same holds to lesser extend for the Th/Ir pair. This 
points to shortcomings in our understanding of the r-process pro- 
duction ratios and/or the chemical history of the star. For CS 
22892-052 all the reference elements except Eu lead to exceed- 
ingly great ages. One might speculate that the Th/Hf age might 
be enlarged by an s-process contribution to the Hf abundance: 
although the star has a very low metalicity, it is also carbon en- 
hanced, so that one cannot rule out the possibility that it may 
have a (so far undetected) companion which went through the 
AGB phase, polluting the photosphere of the primary. However, 
even this hypothesis cannot be invoked to explain the great 
Th/Os and Th/Ir ages. Os and I r are in fact supposed t o show only 
minor ^-process contributions dArlandini et al.l[l999h . Again, we 
must conclude that large systematic uncertainties still exist in the 
theoretical calculations of production ratios f3. 

At this point one may ask what accuracy has been reached in 
spectroscopically measured abundance ratios necessary for ra- 
dioactive dating - in particular for metal-poor stars -, and what 
are the perspectives for improving these accuracies in the future. 
Both questions are not easy to answer. Uncertainties given by 
individual authors are often on the optimistic side. We think that 
a comparison among different authors provides a more robust 
measure of the actually present - as opposed to published - un- 
certainties. Such a comparison might still be biased towards too 
small error bars since different authors often draw from the same 
sources for oscillator strengths, or may use the same observa- 
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Fig. 2. Ages and spectroscopic age uncertainties for star 
CS 22892-052 determined from various chronometer pairs 
(symbols) assuming up to four different pro duction ratio s . Fille d 
symbols refer to the production rations of Kra tz et alj (12 007). 
The d ashed line indicates the age of the universe. Snede n et al.l 
(2003) give a radiochronometric age estimate of 14.2 + 3 Gyr for 
this star. 
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Fig. 3. Same as Fig. [2] for star CS 31082-00 1 In a ddition, for 
this star U/Th ages are available. iHill et ail (120021) give a ra- 
diochronometric age estimate of 14.0 + 2.4 Gyr for this stars. 



tional material. Hence, the analyses are not independent. Sneden 
et al. (l2009h performed such a comparison of the abundances 
which they obtained for the metal- poor r-process en riched halo 
giant CS 3 1082-001 with studies oflHill et all d2002h and Honda 
et al. (l2004h . We repeat this comparison here for the obtained 
abundance ratios. In the comparison Honda-Sneden the list of 
common elements is La, Ce, Pr, Nd, Sm, Eu, Gd, Tb, Dy, Er, 
Tm, Yb, in Hill-Sneden La, Ce, Pr, Nd, Sm, Eu, Gd, Tb, Dy, Er, 
Tm, Hf. Abundance ratios among these elements can be taken as 
a rather good model for the Th/Hf pair since their lines emerge 
from the first, domin ant ionization stage , and from levels of lo w 
excitation potential (lLawler et all 120091: ISneden et ai]l2009h - 
similar to the situation for the Th/Hf pair. Taking the abundance 
ratios from the individual studies ensures that correlation effects 
reducing the errors on abundance ratios are fully taken into ac- 
count. For the rare earth elements often many lines are available 
for the determination of their abundance. This luxury is not af- 
forded by Th or Hf so that we still have a likely bias towards too 
high accuracy when taking the rare earths as example. 

From the lists of elements in common we formed all pos- 
sible abundance ratios (except for an element with itself), and 
calculated the difference of the logarithmic ratios between the 
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various studies. For the pair Hill-Sneden we find a dispersion 
among the differences of the logarithmic ratios of 0.14dex, for 
Honda-Sneden 0.12 dex. In the case of Hill-Sneden, in partic- 
ular the discrepant Hf abundances have to be noted. In the case 
of Honda-Sneden, a substantial part of the dispersion is driven 
by the elements Tm and Tb; leaving them out reduces the dis- 
persion to 0.05 dex; however, this is not suggested by, e.g., ex- 
cessive errors on their abundance as stated by the authors. We 
emphasize that at least in part the various groups use the same 
oscillator strengths. From this exercise we conclude that attribut- 
ing an uncertainty of about 0. 1 dex to the Th/Hf abundance ratio 
is perhaps conservative but not grossly over-estimating the ac- 
tual value. 

How to improve the accuracy? ISneden et al.l (120091) report a 
line-to-line scatter of typically 0.07 dex in a given star when de- 
termining abundances of a rare e arth element and ap plying the 
high-quality oscillator strength of Law ler et al.l ([2009, and refer- 
ences therein). A sub-optimal choice of the model atmosphere 
introduces some scatter, but certainly errors in the oscillator 
strength are an important contributor. Since two abundances are 
involved in the determination of a ratio, an uncertainty of about 
0. 1 dex would result when an individual abundance is given to 
0.07 dex accuracy. It is often argued that by averaging over abun- 
dances of N individual lines one improves the obtained accuracy 
by a factor y/N. However, this would only be the case if the scat- 
ter among the oscillator strength is of random nature which is far 
from obvious. So, even more accurate oscillator strengths of the 
involved line transitions are desirable. This also relates to lines 
of blending species which one has to deal with. At this point an 
appropriate theoretical description of the line formation process 
also comes into play which should account for departures from 
local thermodynamic equilibrium and gas-dynamical effects. All 
this is alleviated by having high signal-to-noise and high resolu- 
tion spectra at hand which ELTs will be able to provide. One 
often has to work with spectra where a line profile is sampled 
with only very few points. Better sampling usually permits more 
accurate measurements. We are convinced that we are not stuck 
at the present level of accuracy to which abundance ratios used in 
radioactive dating can be determined but not insignificant effort 
is necessary to make substantial progress. 



4. The Solar age 

As another application of the dating formulae we tried to deduce 
the solar age. We start by taking the Th/Hf pair as an exam- 
ple since we worke d recently ourselves on this chronometer pair 
(ICaffau et alJ20 08) having all relevant spectroscopic information 
at hand. In the case of the 401.9 nm Thu / 409.3 nm Hfn line 
pair, our favorite spectral synthesis code gives Cr/Cx = 0.13. 
For the Sun one has to consider two complications which are 
less severe for metal-poor stars: 

1 . the Th and Hf observed in the solar system are not the result 
of a single r-process event, but rather of the superposition of 
many events; 

2. while Th is produced only through the r-process, Hf is in 
part produced by the ^-process. 

To take into account the first point one should use an appropri- 
ate model of Galactic chemical evolution (see, e.g., del Peloso et 
al. l2005bh . We ignore this aspect and assume a single r-process 
event, to see if at least a reasonable estimate of the solar age can 
be derived. The second point instead can easily be taken into ac- 
count by using the solar r-process fraction of Hf nuclei. This has 



been estimated bv lLawler et al.1 (120071) to be 56.6 % of the total 
Hf nuclei. For the production ratio we assume ft = 0.862, from 
Table 2 of iKratz et all (120071) . The curve F (^; p = 0.26) related 
to this case is depicted as a dashed line in Fig. Q] For ages 5 and 
12Gyr (i.e. for times corresponding to the age of the Sun and 
to the age of Galactic Globular Clusters, respectively), we find 
aify/t = 4.3 x <t(W r )/W r and 1.9 x cr(W R )/W R , respectively. 
We give these numbers mainly for illustrating the application of 
the formulae involving the equivalent widths. In the solar case 
they are not strictly applicable since when deriving Eq. ( 114l i we 
assumed that uncertainties in the observed spectra are related to 
the photon noise. This is certainly not the case for Solar spec- 
tra. The uncertainties are related to the systematics, in particular 
blending of the thorium line. 

Turning now to spectroscopically measured abundances, 
Table [T] lists the obtained ages and their uncertainties, again for 
different chronometer pairs and assumptions about production 
ratios and r-process fractions (see Tables IA.1I IA.2I IA.3b . The 
Th/Hf age, 22.3 Gyr, comes out to be greater than the age of 
the universe, 13.73 Gyr, as estimat ed from the fluctuat ions of the 
Cosmic Microwave Background dSoergel et alJl2007l). By using 
the A(Hf)=0.88 value suggested bv lLawler et all (l2007h . only 
0.01 dex higher than the abundance we use, the age comes out 
even a bit greater by 0.5 Gyr. Even if the solar abundance of Hf 
and Th were produced in a number of distinct supernova events 
at different times, the resulting mean age as computed from 
Eq. ([8]) cannot be greater than the time elapsed since the ear- 
liest supernova event, which cannot have occurred earlier than 
13.73 Gyr ago. We have ignored the destruction of Th by pho- 
tons in stellar interiors. According to Idel Peloso et akl {2005b) 
the matter cycled through stars may have a reduction of 10-20% 
in Th. However, even when increasing the Solar Th abundance 
by 20%, one still derives an age of 18.2 Gyr from Eq. (0. We 
have to conclude that either the assumed value of ft, or the r- 
process fraction of Hf is incorrect. Among the other chronome- 
ter pairs Th/Eu provides the closest match to the actual Solar 
age, while in general - like for Th/Hf - there it the tendency to 
grossly overestimate the age. The photosphe ric solar abundance 
of europium provided bv lLawler et al.l (120071) is exactly the same 
we are using. Obviously, the complex chemical history of the 
Sun together with uncertainties in the production ratios make it 
very difficult to obtain a reliable Solar age by radioactive dating. 

In this context it is perhaps interesting to note that Sneden 
et al. d2009l) find that the solar photospheric Hf abundance is 
deviating from the meteoritic value in contrast to all its preceed- 
ing (in atomic number) rare earth elements. Similar, the solar 
r-process abundance pattern corresponds to the one of r-process 
enriched iron-poor halo stars except for Hf which appears rela- 
tively depleted in the Sun. Of course, a higher r-process abun- 
dance of Hf in the Sun would even lead to a greater solar age 
aggravating the problem in the radioactive dating. 



5. Conclusions and remarks 

The accuracy of radioactive cosmo-chronometry is not going to 
be extremely high on a single star. A realistic estimate of what 
future instrumentation will be able to achieve does likely not al- 
low one to hope for a better than 2 Gyr precision, taking into 
account observational uncertainties only, leaving aside the sys- 
tematics associated to production rates and chemical enrichment 
models. Although high-mass n-capture elements are currently 
believed to be produced in a fa irly well understood manner by 
the so-called "main" r-process (Farouqi et al. 2008; Kratz et all 
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120071) . the quantitative determination of the production ratios ft 
between the stable and the unstable elements still bear significant 
uncertainties. 

Since in photon noise dominated spectra the measurement 
error on the equivalent width decreases linearly with the observa- 
tional S/N ratio, the attainable dating precision increases roughly 
linearly with the S/N ratio. This constitutes a major practical 
challenge to the application of the method, since many of the 
involved lines are in the blue part of the optical spectrum (e.g. 
for Th, U, and Hf), and the natural target of the analysis are cool 
stars. 

If possible, the decaying element should be chosen such that 
its mean lifetime r is similar to the age one wants to measure. 
This applies also to the "effective" r of a chronometric pair con- 
stituted by two decaying species. From this perspective, the pairs 
U/Th, Th/stable, and U/stable are well suited to measure ages of 
old objects around lOGyr of age. No other naturally occurring 
isotope has a mean lifetime providing a similarly close match. 

Potentially important systematic spectroscopy -related effects 
that we feel deserve further study to fully exploit the potentials 
of radioactive dating are effects related to departures from ther- 
modynamic equilibrium and gas-dynamics on the line formation 
of the chronometric species. For instance, the appli cation of a 
3D rad iation-hydrodynamical model atmosphere bv lCaffau et al. 
(2008) led to a downward revision of the photospheric Solar Th 
abundance by 0. 1 dex. This corresponds to a decrease of Th ages 
by4.7Gyr. 

Globular clusters provide fairly large samples of metal-poor, 
bright, coeval and chemically homogeneous stars. Simultaneous 
measurement of abundances in these stars allows us to reduce 
the observational, statistical uncertainties. This could permit to 
derive r-process production ratios empirically by calibration 
against the known "photometric" age. To our knowledge, this 
has never been attempted so far. 
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Appendix A: Tables of elemental abundances, 
production ratios, and solar r-process fractions 

The following tables list the abundances, production ratios, and 
solar system r-process fractions that we compiled from the liter- 
ature for the purpose of radioactive dating. In Table lATl exclu- 
sively spectroscopically determined (as opposed to meteoritic) 
abundances are give n for the Sun. In T able lA.2l of production ra- 
tios, the values for lKratz et al.l d2007l) refer to the ones listed as 
"Fe-Seed fit2" in Table 2 of their paper. Kratz and collaborators 
report three different estimates of the various production ratios. 



Table A.3. Solar r-process fractions /, for various elements X. 



Table A.2. r-process production ratios. 



Reference 




Production ratio B 






Th/Eu 


Th/Hf Th/Os 


Th/Ir 


U/Th 


Kratz et al. (2007) 


0.453 


0.862 0.093 


0.089 


0.638 


Sneden et al. (2003) 


0.420 








Schatz et al. (2002) 


0.468 


0.071 


0.066 




Cowan et al. (2002) 


0.507 




0.087 





X 


ft 


Reference 


Eu 


0.942 


Arlandini et al. (1999) 


Hf 


0.566 


Lawler et al. (2007) 


Os 


0.92 


Simmerer et al. (2004) 


Ir 


0.99 


Simmerer et al. (2004) 
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Table A.l. r-process elemental abundances in metal-poor stars and the Sun for radioactive dating. 



Object 



Reference 



l^e/HJ 


A(Eu) 




A (del 
A(US) 




A /Thl 

A(in) 


A CT Tl 
A V U ) 


-2.7 


-1.5 








-2.1 




-2.9 


-1.63 ±0.07 




-0.55 ±0.11 




-2.23 ±0.11 


< -2.6 


-3.1 


-0.9 








-1.6 




-3.1 


-0.89 ±0.12 


-0.90 ±0.14 


-0.10 ±0.38 




-1.55 ±0.08 




-3.1 


-0.93 ± 0.09 




-0.05 ±0.10 


+0.00 ± 0.26 


-1.60 ±0.07 




-3.1 


-0.95 ± 0.03 


-0.98 ± 0.10 


+0.02 ±0.13 


-0.10 ±0.10 


-1.57 ± 0.10 


< -2.3 


-2.0 


-0.67 ± 0.05 




+0.45 ±0.10 


+0.30 ±0.13 


— 1.18 ± 0.10 


-2.0 : ±0.30 


-2.18 


-0.86 ± 0.07 


-0.94 ± 0.04 


+0.16 ±0.10 


+0.02 ±0.13 


-1.46 ± 0.05 




-2.90 


-0.63 ± 0.04 


-0.59 ±0.17 


+0.43 ±0.12 


+0.20 ± 0.09 


-0.98 ± 0.02 


-1.92 ±0.11 


-2.32 


-0.88 ± 0.09 








-1.47 ±0.10 




-2.25 


-0.61 ±0.09 








-1.26 ± 0.10 




-2.34 


-1.22 ±0.09 








-1.7 : ±0.10 




-1.23 




+0.08 






-0.68 




-1.29 




-0.02 






-0.63 




-1.22 




+0.03 






-0.83 




-1.30 




-0.12 






-0.88 




-1.20 




+0.08 






-0.48 




-1.25 




+0.03 






-0.68 




-1.23 




-0.02 






-0.58 




-1.22 




+0.03 






-0.63 




-1.09 




+0.18 






-0.68 




-1.28 




-0.12 






-0.58 




-1.33 




-0.22 






-0.68 




-2.38 


-1.48 ±0.02 








-1.99 ±0.07 




-3.15 


-1.81 ±0.03 








-2.36 ± 0.07 




-2.61 


-1.56 ±0.06 








-2.26 ± 0.07 




-2.12 


-1.16 ±0.06 








-1.66 ±0.10 




-2.29 


-1.48 ±0.09 








-1.97 ±0.07 




+0.0 






+ 1.45 ±0.10 


+ 1.35 ±0.10 






+0.0 


+0.52 ± 0.03 












+0.0 




+0.87 ± 0.04 






+0.08 ± 0.03 





HD 1 15444 
HD 1 15444 
CS 22892-052 
CS 22892-052 
CS 22892-052 
CS 22892-052 
BD +17 3248 
HD 221170 
CS 31082-001 
M51 K341 
M51 K462 
M51 K583 
M4L1411 
M4L1501 
M4L1514 
M4 L2406 
M4L2617 
M4 L3209 
M4 L3413 
M4L4511 
M4L4611 
M5 IV-81 
M5 IV-82 
HD 108577 
HD 115444 
HD 186478 
BD +82856 
M92VII-18 
Sun 
Sun 
Sun 



Cowan et al. (1999) 
Westin et al. (2000) 
Cowan et al. (1999) 
Sneden et al. (2000a)" 
Sneden et al. (2000a) 
Sneden et al. (2003) 
Cowan et al. (2002) 
Ivans et al. (2006) 
Hill et al. (2002) 
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a Sneden et al. analysed two different spectra for CS 22892-052 providing two independent measurements. For this reason the object/reference combination is listed twice. 



